Evolutionary remodelling of N-terminal domain loops fine-tunes SARS-CoV-2 spike.

The emergence of SARS-CoV-2 variants has exacerbated the COVID-19 global health crisis. Thus far, all variants carry mutations in the spike glycoprotein, which is a critical determinant of viral transmission being responsible for attachment, receptor engagement and membrane fusion, and an important target of immunity. Variants frequently bear truncations of flexible loops in the N-terminal domain (NTD) of spike; the functional importance of these modifications has remained poorly characterised. We demonstrate that NTD deletions are important for efficient entry by the Alpha and Omicron variants and that this correlates with spike stability. Phylogenetic analysis reveals extensive NTD loop length polymorphisms across the sarbecoviruses, setting an evolutionary precedent for loop remodelling. Guided by these analyses, we demonstrate that variations in NTD loop length, alone, are sufficient to modulate virus entry. We propose that variations in NTD loop length act to fine-tune spike; this may provide a mechanism for SARS-CoV-2 to navigate a complex selection landscape encompassing optimisation of essential functionality, immune-driven antigenic variation and ongoing adaptation to a new host.

Mini-XT, a miniaturized tagmentation-based protocol for efficient sequencing of SARS-CoV-2.

BACKGROUND: The COVID-19 pandemic has highlighted the importance of whole genome sequencing (WGS) of SARS-CoV-2 to inform public health policy. By enabling definition of lineages it facilitates tracking of the global spread of the virus. The evolution of new variants can be monitored and knowledge of specific mutations provides insights into the mechanisms through which the virus increases transmissibility or evades immunity. To date almost 1 million SARS-CoV-2 genomes have been sequenced by members of the COVID-19 Genomics UK (COG-UK) Consortium. To achieve similar feats in a more cost-effective and sustainable manner in future, improved high throughput virus sequencing protocols are required. We have therefore developed a miniaturized library preparation protocol with drastically reduced consumable use and costs. RESULTS: We present the 'Mini-XT' miniaturized tagmentation-based library preparation protocol available on protocols.io ( https://doi.org/10.17504/protocols.io.bvntn5en ). SARS-CoV-2 RNA was amplified using the ARTIC nCov-2019 multiplex RT-PCR protocol and purified using a conventional liquid handling system. Acoustic liquid transfer (Echo 525) was employed to reduce reaction volumes and the number of tips required for a Nextera XT library preparation. Sequencing was performed on an Illumina MiSeq. The final version of Mini-XT has been used to sequence 4384 SARS-CoV-2 samples from N. Ireland with a COG-UK QC pass rate of 97.4%. Sequencing quality was comparable and lineage calling consistent for replicate samples processed with full volume Nextera DNA Flex (333 samples) or using nanopore technology (20 samples). SNP calling between Mini-XT and these technologies was consistent and sequences from replicate samples paired together in maximum likelihood phylogenetic trees. CONCLUSIONS: The Mini-XT protocol maintains sequence quality while reducing library preparation reagent volumes eightfold and halving overall tip usage from sample to sequence to provide concomitant cost savings relative to standard protocols. This will enable more efficient high-throughput sequencing of SARS-CoV-2 isolates and future pathogen WGS.

HDAC6 Degradation Inhibits the Growth of High-Grade Serous Ovarian Cancer Cells.

Histone deacetylase 6 (HDAC6) is a unique histone deacetylating enzyme that resides in the cell cytoplasm and is linked to the modulation of several key cancer related responses, including cell proliferation and migration. The promising anti-cancer response of the first-generation HDAC6 catalytic inhibitors continues to be assessed in clinical trials, although its role in high grade serous ovarian cancer is unclear. This study investigated HDAC6 tumor expression by immunohistochemistry in high-grade serous ovarian cancer (HGSOC) tissue samples and a meta-analysis of HDAC6 gene expression in ovarian cancer from publicly available data. The pharmacological activity of HDAC6 inhibition was assessed in a patient-derived model of HGSOC. HDAC6 was found to be highly expressed in HGSOC tissue samples and in the patient-derived HGSOC cell lines where higher HDAC6 protein and gene expression was associated with a decreased risk of death (hazard ratio (HR) 0.38, (95% confidence interval (CI), 0.16-0.88; p = 0.02); HR = 0.88 (95% CI, 0.78-0.99; p = 0.04)). Similarly, the multivariate analysis of HDAC6 protein expression, adjusting for stage, grade, and cytoreduction/cytoreductive surgery was associated with a decreased risk of death (HR = 0.19 (95% CI, 0.06-0.55); p = 0.002). Knock-down of HDAC6 gene expression with siRNA and protein expression with a HDAC6 targeting protein degrader decreased HGSOC cell proliferation, migration, and viability. Conversely, the selective inhibition of HDAC6 with the catalytic domain inhibitor, Ricolinostat (ACY-1215), inhibited HDAC6 deacetylation of α-tubulin, resulting in a sustained accumulation of acetylated α-tubulin up to 24 h in HGSOC cells, did not produce a robust inhibition of HDAC6 protein function. Inhibition of HGSOC cell proliferation by ACY-1215 was only achieved with significantly higher and non-selective doses of ACY-1215. In summary, we demonstrated, for the first time, that HDAC6 over-expression in HGSOC and all ovarian cancers is a favorable prognostic marker. We provide evidence to suggest that inhibition of HDAC6 catalytic activity with first generation HDAC6 inhibitors has limited efficacy as a monotherapy in HGSOC.

CancerGD: A Resource for Identifying and Interpreting Genetic Dependencies in Cancer.

Genes whose function is selectively essential in the presence of cancer-associated genetic aberrations represent promising targets for the development of precision therapeutics. Here, we present CancerGD, a resource that integrates genotypic profiling with large-scale loss-of-function genetic screens in tumor cell lines to identify such genetic dependencies. CancerGD provides tools for searching, visualizing, and interpreting these genetic dependencies through the integration of functional interaction networks. CancerGD includes different screen types (siRNA, shRNA, CRISPR), and we describe a simple format for submitting new datasets.

RNA-Sequencing data supports the existence of novel VEGFA splicing events but not of VEGFAxxxb isoforms.

Vascular endothelial growth factor (VEGFA), a pivotal regulator of angiogenesis and valuable therapeutic target, is characterised by alternative splicing which generates three principal isoforms, VEGFA121, VEGFA165 and VEGFA189. A second set of anti-angiogenic isoforms termed VEGFAxxxb that utilise an alternative splice site in the final exon have been widely reported, with mRNA detection based principally upon RT-PCR assays. We sought confirmation of the existence of the VEGFAxxxb isoforms within the abundant RNA sequencing data available publicly. Whilst sequences derived specifically from each of the canonical VEGFA isoforms were present in many tissues, there were no sequences derived from VEGFAxxxb isoforms. Sequencing of approximately 50,000 RT-PCR products spanning the exon 7-8 junction in 10 tissues did not identify any VEGFAxxxb transcripts. The absence or extremely low expression of these transcripts in vivo indicates that VEGFAxxxb isoforms are unlikely to play a role in normal physiology. Our analyses also revealed multiple novel splicing events supported by more reads than previously reported for VEGFA145 and VEGFA148 isoforms, including three from novel first exons consistent with existing transcription start site data. These novel VEGFA isoforms may play significant roles in specific cell types.

Sequence and Expression of Complement Factor H Gene Cluster Variants and Their Roles in Age-Related Macular Degeneration Risk.

PURPOSE: To investigate how potentially functional genetic variants are coinherited on each of four common complement factor H (CFH) and CFH-related gene haplotypes and to measure expression of these genes in eye and liver tissues. METHODS: We sequenced the CFH region in four individuals (one homozygote for each of four common CFH region haplotypes) to identify all genetic variants. We studied associations between the haplotypes and AMD phenotypes in 2157 cases and 1150 controls. We examined RNA-seq profiles in macular and peripheral retina and retinal pigment epithelium/choroid/sclera (RCS) from eight eye donors and three liver samples. RESULTS: The haplotypic coinheritance of potentially functional variants (including missense variants, novel splice sites, and the CFHR3-CFHR1 deletion) was described for the four common haplotypes. Expression of the short and long CFH transcripts differed markedly between the retina and liver. We found no expression of any of the five CFH-related genes in the retina or RCS, in contrast to the liver, which is the main source of the circulating proteins. CONCLUSIONS: We identified all genetic variants on common CFH region haplotypes and described their coinheritance. Understanding their functional effects will be key to developing and stratifying AMD therapies. The small scale of our expression study prevented us from investigating the relationships between CFH region haplotypes and their expression, and it will take time and collaboration to develop epidemiologic-scale studies. However, the striking difference between systemic and ocular expression of complement regulators shown in this study suggests important implications for the development of intraocular and systemic treatments.

Novel insights into the insect trancriptome response to a natural DNA virus.

BACKGROUND: Little is known about invertebrate responses to DNA viruses. Here, we infect a commercially important pest moth species Plodia interpunctella with its naturally infecting DNA virus. We sequenced, assembled and annotated the complete transcriptome of the moth, and a partial transcriptome of the virus. We then tested for differential gene expression between moths that were exposed to the virus and controls. RESULTS: We found 51 genes that were differentially expressed in moths exposed to a DNA baculovirus compared to controls. Gene set enrichment analysis revealed that cuticle proteins were significantly overrepresented in this group of genes. Interestingly, 6 of the 7 differentially expressed cuticle proteins were downregulated, suggesting that baculoviruses are able to manipulate its host's response. In fact, an additional 29 of the 51 genes were also downregulated in exposed compared with control animals, including a gram-negative binding protein. In contrast, genes involved in transposable element movement were upregulated after infection. CONCLUSIONS: We present the first experiment to measure genome-wide gene expression in an insect after infection with a natural DNA virus. Our results indicate that cuticle proteins might be key genes underpinning the response to DNA viruses. Furthermore, the large proportion of genes that were downregulated after viral exposure suggests that this virus is actively manipulating the insect immune response. Finally, it appears that transposable element activity might increase during viral invasion. Combined, these results provide much needed host candidate genes that respond to DNA viral invaders.

Small RNAs from plants, bacteria and fungi within the order Hypocreales are ubiquitous in human plasma.

BACKGROUND: The human microbiome plays a significant role in maintaining normal physiology. Changes in its composition have been associated with bowel disease, metabolic disorders and atherosclerosis. Sequences of microbial origin have been observed within small RNA sequencing data obtained from blood samples. The aim of this study was to characterise the microbiome from which these sequences are derived. RESULTS: Abundant non-human small RNA sequences were identified in plasma and plasma exosomal samples. Assembly of these short sequences into longer contigs was the pivotal novel step in ascertaining their origin by BLAST searches. Most reads mapped to rRNA sequences. The taxonomic profiles of the microbes detected were very consistent between individuals but distinct from microbiomes reported at other sites. The majority of bacterial reads were from the phylum Proteobacteria, whilst for 5 of 6 individuals over 90% of the more abundant fungal reads were from the phylum Ascomycota; of these over 90% were from the order Hypocreales. Many contigs were from plants, presumably of dietary origin. In addition, extremely abundant small RNAs derived from human Y RNAs were detected. CONCLUSIONS: A characteristic profile of a subset of the human microbiome can be obtained by sequencing small RNAs present in the blood. The source and functions of these molecules remain to be determined, but the specific profiles are likely to reflect health status. The potential to provide biomarkers of diet and for the diagnosis and prognosis of human disease is immense.

Quality control of next-generation sequencing data without a reference.

Next-generation sequencing (NGS) technologies have dramatically expanded the breadth of genomics. Genome-scale data, once restricted to a small number of biomedical model organisms, can now be generated for virtually any species at remarkable speed and low cost. Yet non-model organisms often lack a suitable reference to map sequence reads against, making alignment-based quality control (QC) of NGS data more challenging than cases where a well-assembled genome is already available. Here we show that by generating a rapid, non-optimized draft assembly of raw reads, it is possible to obtain reliable and informative QC metrics, thus removing the need for a high quality reference. We use benchmark datasets generated from control samples across a range of genome sizes to illustrate that QC inferences made using draft assemblies are broadly equivalent to those made using a well-established reference, and describe QC tools routinely used in our production facility to assess the quality of NGS data from non-model organisms.

Biomineral proteins from Mytilus edulis mantle tissue transcriptome.

The common blue mussel, Mytilus edulis, has a bimineralic shell composed of approximately equal proportions of the two major polymorphs of calcium carbonate: calcite and aragonite. The exquisite biological control of polymorph production is the focus of research interest in terms of understanding the details of biomineralisation and the proteins involved in the process of complex shell formation. Recent advances in ease and availability of pyrosequencing and assembly have resulted in a sharp increase in transcriptome data for invertebrate biominerals. We have applied Roche 454 pyrosequencing technology to profile the transcriptome for the mantle tissue of the bivalve M. edulis. A comparison was made between the results of several assembly programs: Roche Newbler assembler versions 2.3, 2.5.2 and 2.6 and MIRA 3.2.1 and 3.4.0. The Newbler and MIRA assemblies were subsequently merged using the CAP3 assembler to give a higher consensus in alignments and a more accurate estimate of the true size of the M. edulis transcriptome. Comparison sequence searches show that the mantle transcripts for M. edulis encode putative proteins exhibiting sequence similarities with previously characterised shell proteins of other species of Mytilus, the Bivalvia Pinctada and haliotid gastropods. Importantly, this enhanced transcriptome has detected several transcripts that encode proteins with sequence similarity with previously described shell biomineral proteins including Shematrins and lysine-rich matrix proteins (KRMPs) not previously found in Mytilus.

Secretion of protective antigens by tissue-stage nematode larvae revealed by proteomic analysis and vaccination-induced sterile immunity.

Gastrointestinal nematode parasites infect over 1 billion humans, with little evidence for generation of sterilising immunity. These helminths are highly adapted to their mammalian host, following a developmental program through successive niches, while effectively down-modulating host immune responsiveness. Larvae of Heligmosomoides polygyrus, for example, encyst in the intestinal submucosa, before emerging as adult worms into the duodenal lumen. Adults release immunomodulatory excretory-secretory (ES) products, but mice immunised with adult H. polygyrus ES become fully immune to challenge infection. ES products of the intestinal wall 4th stage (L4) larvae are similarly important in host-parasite interactions, as they readily generate sterile immunity against infection, while released material from the egg stage is ineffective. Proteomic analyses of L4 ES identifies protective antigen targets as well as potential tissue-phase immunomodulatory molecules, using as comparators the adult ES proteome and a profile of H. polygyrus egg-released material. While 135 proteins are shared between L4 and adult ES, 72 are L4 ES-specific; L4-specific proteins correspond to those whose transcription is restricted to larval stages, while shared proteins are generally transcribed by all life cycle forms. Two protein families are more heavily represented in the L4 secretome, the Sushi domain, associated with complement regulation, and the ShK/SXC domain related to a toxin interfering with T cell signalling. Both adult and L4 ES contain extensive but distinct arrays of Venom allergen/Ancylostoma secreted protein-Like (VAL) members, with acetylcholinesterases (ACEs) and apyrase APY-3 particularly abundant in L4 ES. Serum antibodies from mice vaccinated with L4 and adult ES react strongly to the VAL-1 protein and to ACE-1, indicating that these two antigens represent major vaccine targets for this intestinal nematode. We have thus defined an extensive and novel repertoire of H. polygyrus proteins closely implicated in immune modulation and protective immunity.

The transcriptome of the invasive eel swimbladder nematode parasite Anguillicola crassus.

BACKGROUND: Anguillicola crassus is an economically and ecologically important parasitic nematode of eels. The native range of A. crassus is in East Asia, where it infects Anguilla japonica, the Japanese eel. A. crassus was introduced into European eels, Anguilla anguilla, 30 years ago. The parasite is more pathogenic in its new host than in its native one, and is thought to threaten the endangered An. anguilla across its range. The molecular bases for the increased pathogenicity of the nematodes in their new hosts is not known. RESULTS: A reference transcriptome was assembled for A. crassus from Roche 454 pyrosequencing data. Raw reads (756,363 total) from nematodes from An. japonica and An. anguilla hosts were filtered for likely host contaminants and ribosomal RNAs. The remaining 353,055 reads were assembled into 11,372 contigs of a high confidence assembly (spanning 6.6 Mb) and an additional 21,153 singletons and contigs of a lower confidence assembly (spanning an additional 6.2 Mb). Roughly 55% of the high confidence assembly contigs were annotated with domain- or protein sequence similarity derived functional information. Sequences conserved only in nematodes, or unique to A. crassus were more likely to have secretory signal peptides. Thousands of high quality single nucleotide polymorphisms were identified, and coding polymorphism was correlated with differential expression between individual nematodes. Transcripts identified as being under positive selection were enriched in peptidases. Enzymes involved in energy metabolism were enriched in the set of genes differentially expressed between European and Asian A. crassus. CONCLUSIONS: The reference transcriptome of A. crassus is of high quality, and will serve as a basis for future work on the invasion biology of this important parasite. The polymorphisms identified will provide a key tool set for analysis of population structure and identification of genes likely to be involved in increased pathogenicity in European eel hosts. The identification of peptidases under positive selection is a first step in this programme.

Transcriptional changes in Teladorsagia circumcincta upon encountering host tissue of differing immune status.

The aim of this study was to elucidate transcriptional changes in the parasitic nematode Teladorsagia circumcincta upon encountering either naïve or immune ovine hosts. Pools of 100 000 exsheathed 3rd- stage T. circumcincta larvae were exposed in vitro to either an immune or naïve ovine abomasal environment, RNA was extracted from the larvae and sequenced using the Roche 454 platform. Each sample produced approximately 82 000 reads that assembled to give approximately 5500 Isotigs (contigs). The two sequence datasets were clustered together to give a total of 6969 clusters of which 18 were differentially expressed (P<0·001) between the two groups. Clusters with a predominance of reads in larvae exposed to the immune abomasal environment encoded homologues of peptidyl-glycine alpha-amidating monooxygenase, heat shock-protein 16-2 and IDA-1, a tyrosine phosphatase-like receptor protein. Clusters with a predominance of reads in the naïve environment encoded homologues of cytochrome b, EGg Laying defective family member 21 and NADH dehydrogenase subunit 5. Gene ontology analyses indicated that larvae exposed to the immune environment showed an increase in expression of genes involved in 'carbon utilization', 'response to stimulus' and 'developmental process'. These data suggest that T. circumcincta modulates gene expression in response to the immune status of the host.

Polyphenism in social insects: insights from a transcriptome-wide analysis of gene expression in the life stages of the key pollinator, Bombus terrestris.

BACKGROUND: Understanding polyphenism, the ability of a single genome to express multiple morphologically and behaviourally distinct phenotypes, is an important goal for evolutionary and developmental biology. Polyphenism has been key to the evolution of the Hymenoptera, and particularly the social Hymenoptera where the genome of a single species regulates distinct larval stages, sexual dimorphism and physical castes within the female sex. Transcriptomic analyses of social Hymenoptera will therefore provide unique insights into how changes in gene expression underlie such complexity. Here we describe gene expression in individual specimens of the pre-adult stages, sexes and castes of the key pollinator, the buff-tailed bumblebee Bombus terrestris. RESULTS: cDNA was prepared from mRNA from five life cycle stages (one larva, one pupa, one male, one gyne and two workers) and a total of 1,610,742 expressed sequence tags (ESTs) were generated using Roche 454 technology, substantially increasing the sequence data available for this important species. Overlapping ESTs were assembled into 36,354 B. terrestris putative transcripts, and functionally annotated. A preliminary assessment of differences in gene expression across non-replicated specimens from the pre-adult stages, castes and sexes was performed using R-STAT analysis. Individual samples from the life cycle stages of the bumblebee differed in the expression of a wide array of genes, including genes involved in amino acid storage, metabolism, immunity and olfaction. CONCLUSIONS: Detailed analyses of immune and olfaction gene expression across phenotypes demonstrated how transcriptomic analyses can inform our understanding of processes central to the biology of B. terrestris and the social Hymenoptera in general. For example, examination of immunity-related genes identified high conservation of important immunity pathway components across individual specimens from the life cycle stages while olfactory-related genes exhibited differential expression with a wider repertoire of gene expression within adults, especially sexuals, in comparison to immature stages. As there is an absence of replication across the samples, the results of this study are preliminary but provide a number of candidate genes which may be related to distinct phenotypic stage expression. This comprehensive transcriptome catalogue will provide an important gene discovery resource for directed programmes in ecology, evolution and conservation of a key pollinator.

Gene conversion rapidly generates major histocompatibility complex diversity in recently founded bird populations.

Population bottlenecks can restrict variation at functional genes, reducing the ability of populations to adapt to new and changing environments. Understanding how populations generate adaptive genetic variation following bottlenecks is therefore central to evolutionary biology. Genes of the major histocompatibility complex (MHC) are ideal models for studying adaptive genetic variation due to their central role in pathogen recognition. While de novo MHC sequence variation is generated by point mutation, gene conversion can generate new haplotypes by transferring sections of DNA within and across duplicated MHC loci. However, the extent to which gene conversion generates new MHC haplotypes in wild populations is poorly understood. We developed a 454 sequencing protocol to screen MHC class I exon 3 variation across all 13 island populations of Berthelot's pipit (Anthus berthelotii). We reveal that just 11-15 MHC haplotypes were retained when the Berthelot's pipit dispersed across its island range in the North Atlantic ca. 75,000 years ago. Since then, at least 26 new haplotypes have been generated in situ across populations. We show that most of these haplotypes were generated by gene conversion across divergent lineages, and that the rate of gene conversion exceeded that of point mutation by an order of magnitude. Gene conversion resulted in significantly more changes at nucleotide sites directly involved with pathogen recognition, indicating selection for functional variants. We suggest that the creation of new variants by gene conversion is the predominant mechanism generating MHC variation in genetically depauperate populations, thus allowing them to respond to pathogenic challenges.

Proteomic analysis of secretory products from the model gastrointestinal nematode Heligmosomoides polygyrus reveals dominance of venom allergen-like (VAL) proteins.

The intestinal helminth parasite, Heligmosomoides polygyrus bakeri offers a tractable experimental model for human hookworm infections such as Ancylostoma duodenale and veterinary parasites such as Haemonchus contortus. Parasite excretory-secretory (ES) products represent the major focus for immunological and biochemical analyses, and contain immunomodulatory molecules responsible for nematode immune evasion. In a proteomic analysis of adult H. polygyrus secretions (termed HES) matched to an extensive transcriptomic dataset, we identified 374 HES proteins by LC-MS/MS, which were distinct from those in somatic extract HEx, comprising 446 identified proteins, confirming selective export of ES proteins. The predominant secreted protein families were proteases (astacins and other metalloproteases, aspartic, cysteine and serine-type proteases), lysozymes, apyrases and acetylcholinesterases. The most abundant products were members of the highly divergent venom allergen-like (VAL) family, related to Ancylostoma secreted protein (ASP); 25 homologues were identified, with VAL-1 and -2 also shown to be associated with the parasite surface. The dominance of VAL proteins is similar to profiles reported for Ancylostoma and Haemonchus ES products. Overall, this study shows that the secretions of H. polygyrus closely parallel those of clinically important GI nematodes, confirming the value of this parasite as a model of helminth infection.

Comparative genome analysis provides insights into the evolution and adaptation of Pseudomonas syringae pv. aesculi on Aesculus hippocastanum.

A recently emerging bleeding canker disease, caused by Pseudomonas syringae pathovar aesculi (Pae), is threatening European horse chestnut in northwest Europe. Very little is known about the origin and biology of this new disease. We used the nucleotide sequences of seven commonly used marker genes to investigate the phylogeny of three strains isolated recently from bleeding stem cankers on European horse chestnut in Britain (E-Pae). On the basis of these sequences alone, the E-Pae strains were identical to the Pae type-strain (I-Pae), isolated from leaf spots on Indian horse chestnut in India in 1969. The phylogenetic analyses also showed that Pae belongs to a distinct clade of P. syringae pathovars adapted to woody hosts. We generated genome-wide Illumina sequence data from the three E-Pae strains and one strain of I-Pae. Comparative genomic analyses revealed pathovar-specific genomic regions in Pae potentially implicated in virulence on a tree host, including genes for the catabolism of plant-derived aromatic compounds and enterobactin synthesis. Several gene clusters displayed intra-pathovar variation, including those encoding type IV secretion, a novel fatty acid biosynthesis pathway and a sucrose uptake pathway. Rates of single nucleotide polymorphisms in the four Pae genomes indicate that the three E-Pae strains diverged from each other much more recently than they diverged from I-Pae. The very low genetic diversity among the three geographically distinct E-Pae strains suggests that they originate from a single, recent introduction into Britain, thus highlighting the serious environmental risks posed by the spread of an exotic plant pathogenic bacterium to a new geographic location. The genomic regions in Pae that are absent from other P. syringae pathovars that infect herbaceous hosts may represent candidate genetic adaptations to infection of the woody parts of the tree.

InCa-SiteFinder: a method for structure-based prediction of inositol and carbohydrate binding sites on proteins.

Carbohydrate binding sites are considered important for cellular recognition and adhesion and are important targets for drug design. In this paper we present a new method called InCa-SiteFinder for predicting non-covalent inositol and carbohydrate binding sites on the surface of protein structures. It uses the van der Waals energy of a protein-probe interaction and amino acid propensities to locate and predict carbohydrate binding sites. The protein surface is searched for continuous volume envelopes that correspond to a favorable protein-probe interaction. These volumes are subsequently analyzed to demarcate regions of high cumulative propensity for binding a carbohydrate moiety based on calculated amino acid propensity scores. InCa-SiteFinder(1) was tested on an independent test set of 80 protein-ligand complexes. It efficiently identifies carbohydrate binding sites with high specificity and sensitivity. It was also tested on a second test set of 80 protein-ligand complexes containing 40 known carbohydrate binders (having 40 carbohydrate binding sites) and 40 known drug-like compound binders (having 58 known drug-like compound binding sites) for the prediction of the location of the carbohydrate binding sites and to distinguish these from the drug-like compound binding sites. At 73% sensitivity the method showed 98% specificity. Almost all of the carbohydrate and drug-like compound binding sites were correctly identified with an overall error rate of 12%.